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Abstract 



We investigate several spin-related phenomena in a confined two-dimensional 



electron gas (2DEG) using the Hartree-Fock approximation for the mutual 

a 

i . Coulomb interaction of the electrons. The exchange term of the interaction 

causes a large splitting of the spin levels whenever the chemical potential 
lies within a Landau band (LB). This splitting can be reinterpreted as an 



^ . enhancement of an effective g- factor, g* . The increase of g* when a LB is 

half filled can be accompanied by a spontaneous formation of a static spin- 
inversion state (SIS) whose details depend on the system size. The coupling 
of the states of higher LB's into the lowest band by the Coulomb interaction 
of the 2DEG is essential for the SIS to occur. The far-infrared absorption of 
the system, relatively insensitive to the spin splitting, develops clear signs of 
the SIS. 
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I. INTRODUCTION 



The effects of the exchange interaction on the appearance of macroscopic 
spin structures have been studied in semiconductor microstructures in re- 
duced dimensions by several researchers both theoretically and in experi- 
ments. The enhancement of the effective g-factor, g* , of a two-dimensional 
electron gas (2DEG) in the quantum Hall regime has been reviewed by Ando, 
Fowler, and Stern.0 For the unbounded 2DEG Ando and Uemurai presented 
a model where the broadening of the Landau levels due to impurity scattering 
is treated in the self-consistent Born approximation (SCBA). The dielectric 
function is calculated with the inclusion of the lowest order exchange energy 
of the screened Coulomb interaction in the self-energy of the electrons. For a 
strong magnetic field the overlapping of Landau levels with different indices 
is neglected. 

The enhancement of g* can lead to a spin polarization of the 2DEG at 
certain values of the filling factor v, and, in addition, the exchange interac- 
tion can lead to the spontaneous formation of spin-density! or charge-density 
wavesBi. The onset of a spin-density wave state in a parabolic quantum well 
has been studied by Brey and Halperin using a modified Hartree-Fock ap- 
proximation (HFA) with a point-contact exchange interaction. They find a 
divergence of the electric susceptibility in the presence of a magnetic field of 
intermediate strength parallel to the quantum well and an infinitesimal fic- 
titous magnetic field perpendicular to the quantum wellJl This spin-density 
wave state has a wavevector along the quantum well parallel to the interme- 
diate magnetic field and occurs only when the quantum well is wide enough 
and the exchange interaction has a strength larger than a critical value. The 
calculated optical properties of a <5-doped quantum well in the HFA due to 
spin- and charge-density excitations have been found to be in good agree- 
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ment with experiments,! as well as those of donor states in 2DEG in strong 
magnetic fields.! In the quantum well Hembree et aZ.i discovered abrupt spin 
polarization of the system at high magnetic fields and a spin-inversion regime 
where the net spin alignment strongly varies across the well. They studied 
the effect in different approximation schemes and in the presence of impu- 
rity scattering. Recently the effects of the g-factor enhancement on various 
transport coefficients has been reported by the same group.0 

As to microstructures of further reduced dimensionality the spontaneous 
polarization of of an array of quantum dots into a ferroelectric or antiferro- 
electric state has been investigated by Kempa, Broido, and Bakshi.0 The 
spin degree of freedom together with the exchange interaction and correlation 
effects have also been found to be essential to model few electrons in a single 
quantum dot in magnetic field.BHll 

In this paper we are concerned with the spin-related phenomena associated 
with the exchange interaction that can occur in quantum dots with a large 
number of electrons. We study the spin splitting of Landau bands (LB's) due 
to the enhancement of g* , and the formation of a spin-inversion state (SIS) 
in a strictly two-dimensional finite size electron system in a perpendicular 
magnetic field of intermediate strength. The system size is chosen to be of 



the order of several magnetic lengths, I = y/hc/(eB). The LB's in the center 
of the system do approach flat Landau levels indicating that an electron in 
the center does not feel the boundary. We are thus able to study the crossing 
from the quantum regime in which the electronic confinement dominates over 
the electron-electron interaction to the regime in which electrostatics plays a 
dominant role. Finally we show how the formation of a SIS can be detected 
in the far-infrared (FIR) absorption spectrum of the system. 
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II. MODEL 



We consider N s strictly two-dimensional electrons to model qualitatively 
a real heterostructure where the 2DEG is confined to the lowest electrical 
subband. The 2DEG is confined to a disk of radius R in the 2D-plane by a 
potential step 



V conf( r ) = 



exp 



R-r 
4Ar 



+ 1 



where Ar = 22 A. To ensure charge neutrality of the system a positive back- 
ground charge +enb resides on the disk 



n h (r) = n s 



exp 



r-R 
Ar 



+ 1 



with the average electron density of the system given by n s = N s /(irR 2 ) = 
(n s (r)). In the HFA the state of each electron is described by a single-electron 
Schrodinger equation 

{H + V H (r) + V con t(r)}Mr) 
-J d 2 r' A(f,f')Mr') = e a M^) 

for an electron moving in a Hartree potential 

- '■ ,n s (r') - n b (r') 



V H (r) = — f d 2 r' 

K J 



caused by the charge density —e{n s (r) — n^r)}, and a nonlocal exchange 
potential with 



A(f,f') = -£/(^-^) 



r — r ' 



The equilibrium occupation of the electronic states is according to the Fermi 
distribution f{ep — fi) at finite temperature T. The density of the electrons 
n s (r) is constructed from the energy spectrum {e a } and the wave functions 



,(r) = El^a(*Ol a /(^-M), 



together with the chemical potential [i. The label a represents the radial 
quantum number n r , the angular quantum number M, and the spin quantum 
number s = H° is the single particle Hamiltonian for one electron with 
spin in a constant perpendicular external magnetic field JlZHli A Landau band 
index n can be constructed from the quantum numbers n r and M as n = 
(\M\ — M)/2 + n r . The Landau levels of H° with energy E nj M,s = fioj c {n + 
\) + sg*(nB/ft)B are degenerate with respect to M with the degeneracy no = 
(2-7r/ 2 ) _1 per spin orientation. fiB is the Bohr magneton (eh/2mc). The 
cyclotron frequency is given by uj c = eB/(mc). The Hartree-Fock energy 
spectrum {e a } and the corresponding wave functions {V'q} ar e now found by 
solving ©-(H) iteratively in the basis of H°, The chemical potential \i is 

recalculated in each iteration in order to preserve the total number of electrons 
N s . The number of basis functions used in the diagonalization is chosen 
such that a further increase of the subset results in an unchanged density 
n s (r). The calculations have been repeated for several initial conditions with 
different spin configuration in order to search for the ground state of the 
Hartree-Fock equations The total energy of the system E^ Q ^ can be 

found by summing up the single electron contributions and carefully counting 
the interaction energy of each electron pair only onceJll 

The FIR-absorption of the system is calculated as a self-consistent linear 
response to an external potential^ 

(f xt (f, t) = £ ext r exp {-iNpip - i{oo + irf)t} , 

where rj — > + . N p = ±1 corresponds to left or right circular polarization. 
The small size of the system compared to the wavelength of the external 
radiation makes possible to use a electrostatic potential representing a time 
dependent but spatially constant external electrical field ^ ex ^ = — V0 ex ^. 



In this so-called time-dependent HFA the change of the density matrix due 
to an adiabatically switched-on total electrostatic potential </> sc is calculated 
within a linear approximation. The total potential consists of the external 
potential and the induced potential cj)^ n( ^ = <ft H + (fi F due to the direct and the 
exchange interaction of the electrons. The induced potential in turn depends 
on the density matrix, thus closing the circle and allowing for a self-consistent 
evaluation of the total potential together with an expression for the frequency 
dependent dielectric tensor e a p^(u). The power absorption is then calculated 
from the Joule heating of the system due to ^ ex ^ 



P(u) = e£ ext J2 {El3 h jE " ) (/?|ria)27rJ M/3 ,Af Q±j v^{r /3 (^)(a|(-e ( A SC )|/3)} , 

a/3 

where £ ex ^ is the strength of the external field and 

f#(u) 



ea/3/, Vi _ 1 J fy f° 



H }oj + (up — u a ) + irj 
with the Fermi distribution f a = f(e a — fj). 



III. RESULTS 

The calculations for the box-like confinement ([l]) are carried out with 
GaAs parameters: m* = 0.067m e , n = 12.4, and g* = —0.44. The occupation 
of the LB's is varied by changing the number of electrons N s at a constant 
strength of the magnetic field B = 3.0 T (we could equally have changed the 
magnetic field keeping constant the number of particles). Since the radius 
of the system R > 1000 A is much larger than the magnetic length I ~ 
148 A and the effective Bohr radius ag ~ 97.9 A we can use as a convenient 
label an effective filling factor v describing the occupation of the LB's in the 
interior of the system. The cyclotron energy hu c ~ 5.2 meV, so a sufficient 
height of the confining potential is Uq = 60 meV in order to include several 
LB's in the calculation. For B = 3.0 T the bare spin splitting of the LB's 
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-^Zeeman = {.9* I^B/fyB ~ 0.076 meV is much smaller than their separation 
huj c and corresponds to the thermal energy k^T at T ~ 0.9 K. 

Figure [l] shows the HFA quasiparticle energy spectrum for four values of 
N s such that the filling factor v (defined as the number of occupied bands in 
the central region of the box) ranges from 4 to 2. Figure [l]a corresponds to 
the case v = 4. Electrons in the first LB (n = 0) form a large paramagnetic 
compact droplet while electrons in the second LB (n = 1) form a smaller 
one. Figures |]b and [pc show clearly a large spin splitting of the LB's due 
to the enhancement of g* when the 2nd LB is half filled {y = 3) and the 
electrons in it form a ferromagnetic compact droplet. Finally, Fig. ||d shows 
the case corresponding to the droplet at v = 2 when no electrons are left in 
the 2nd LB. In our case the 2nd LB behaves in all respects like an indepen- 
dent, smaller quantum dot, and its properties are identical to those studied 
previously.il! However, the first LB present a more complicated behavior: The 
energy spectra for six values of N s such that the chemical potential lies in 
the neighborhood of the first LB is seen in Fig. Now the filling factor lies 
within the range 1 < v < 2. For N s = 42 (Fig. |2|a) both spin bands of the 
first LB in the bulk region are still filled (fj, is still lying between the first 
and second LB's but closer to the former one). The small Zeeman energy 
makes the LB's look degenerate with respect to the spin degree of freedom. 
As the number of electrons is r educed to 38 (Fig. |b) the spin bands split 
up near the edge and the number of spin-down electrons becomes smaller 
than that of spin-up electrons. In addition to this splitting (which was also 
present in the 2nd LB in Fig. [I]) we can observe two instability points (one 
for each spin band) rising near the center of the system. By instability points 
we mean bumps in the spin bands approaching the chemical potential. Thus, 
one should expect some transition as the number of electrons keeps changing 
(or the magnetic field) and these bumps touch the chemical potential. Before 

7 



that can happen one can even see signatures of such transitions in the density 
plotted in Fig. [|. The finite temperature "reveals" a budding spin-inversion 
state due to the difference in distance to the chemical potential between the 
spin- up and spin-down bands for a given position (M). The onset of such 
a SIS takes place when the number of particles is reduced further and those 
bumps cross the chemical potential (Figs. §c to |2|e). Finally, the compact 
droplet at v = 1 is formed (Fig. ||f). 

In order to understand better the formation of the SIS let us consider 
a simpler but equivalent situation. The number of LB's cannot be reduced 
considerably in the calculations with the box-like confinement (||) since we are 
using a basis constructed of the eigenfunctions of noninteracting electrons in 
an infinite system. However, by considering parabolically confined interacting 
electrons and using the one-electron basis set of such a system we are left 
only with band mixing due to the electron-electron interaction. We have thus 
calculated the energy spectra of a 2DEG in a parabolic quantum dot with 
confinement frequency Uujq for an increasing number of LB's (from one to 
three) at T = 1.0K. First, we analyze the case of one LB at zero temperature. 
Figure || shows the evolution with B of the band structure for N s = 30. The 
spin splitting opens up from the edge to the center of the LB in the parabolic 
confinement as we go from v = 2 to v = 1. Surprisingly, one can see how both 
spin bands near the center of the system bend upward, and, eventually, one of 
them crosses the chemical potential. This cannot happen for a smaller number 
of electrons since in that case the s = — \ electrons leave their band before 
the unstable point (in the center) crosses the chemical potential. Such central 
instability requires a certain size of the electronic droplet and constitutes the 
initial stage of the SIS. If we include higher LB's the spectra becomes more 
complicated and the instability points of each spin band shift from each other 
due to the mixing with higher LB's. This result is presented in Fig. [|. The 
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total energy of the system -E^ot * s 1789.07 meV for one LB, 1728.3 meV for 
two LB's, and 1718.5 meV for three. Obviously the calculation with one LB 
does not represent well the ground state for the given values of hco c and hujQ, 
but helps us to get an insight on the spin instability. No twisting of the spin 
bands is ever seen for the calculation with one LB. This can be verified by 
checking the analytical expressions for the matrix elements of the exchange 
interactions. The additional degree of freedom introduced to the system by 
allowing coupling of states of higher LB's into the lowest LB for interacting 
electrons is essential in order to obtain the full richness of the spin band 
structure. 

The formation of the SIS invokes clear signs in the FIR spectrum P(to) of 
the 2DEG detailed in Fig. |§ The first two subfi gures show the spectrum in 
the Hartree approximation (HA) and the HFA, respectively. In the HA the ex- 
change interaction is neglected both in the ground state and the excited states. 
A common feature is the occurrence of two strong absorption lines, the lower 
one in energy corresponding to N p = +1 and the higher one corresponding to 
N p = — 1. These two lines can either be identified as the ones corresponding 
to the center of mass motion predicted by the generalized Kohn theorem for 
quantum dots with parabolic confinement or more appropriately here as 
the low energy excitation of an edge plasmon and the 2D bulk plasmon at 
energy slightly higher than the cyclotron resonance E c = TlujJ^ Both approx- 
imation then show small absorption peaks above the bulk magnetoplasmon 
that have been identified as absorption due to single electron transitions .00 
The spin splitting itself does not have large effects on the absorption due to 
the bulk magnetoplasmon but the finer details of the corresponding absorp- 
tion peak in a parabolic quantum well have been studied by Hembree et. al.,0 
here we shall concentrate on the effects of the SIS. By comparing the spectra 
for the two approximations at energy below the energy of the edge plasmon 

9 



we find small peaks for N p = — 1 that are enlargened in the last subfigure of 
Fig. ||. No such peaks are found in the HA. They are only present when the 
SIS occurs and the ones with the lowest energy are caused by single electron 
transitions in the lowest LB, intra-Landau-band transitions with M — > M — 1 
that are only possible because of the twisting of the LB's. Corresponding 
absorption peaks of the opposite polarization N p = +1 can also be found in 
the center subfigure at similar energy, but the peaks with N p = — 1 are much 
more characteristic of the SIS since otherwise peaks of that polarization are 
never found for low energy. As soon as the spin Landau bands of the lowest 
LB cross twice a second absorption peak appears with energy above the edge 
magnetoplasmon but below the bulk plasmon. The occurrence of this second 
row of peaks has to be correlated with the fact that the twisting of the lowest 
and the next LB, that did mirror one another for lower N s , are now out of 
phase for the higher values of N s corresponding to v just below 2. 



IV. DISCUSSION AND SUMMARY 

In a system of a confined 2DEG we have been able to demonstrate both 
bulk effects and phenomena caused by the finite size of the system, in the 
absence of any impurity scattering of the electrons. The 2D system is large 
enough so that the LB's approach flat Landau levels for low values of the 
angular quantum number M. This can be interpreted as the formation of 2D 
bulk states inside the system. The ensuing singular density of states together 
with the exchange interaction causes the well known oscillations of the energy 
separation of the LB's with the same Landau level index n but opposite spin 
orientations as a function of the filling factor v. Here we have seen that the 
enhancement of g* occurs not only in the LB where \x is located but in all the 
LB's included in the model. Similar behavior has been established in optical 
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measurements of a 2DEG by Kukushkin.EH 

We have observed the spontaneous formation of concentric circular regions 
of different spin phases when the spin splitting of the first LB's is opening up 
with a decreasing v at a low temperature. The shape of this SIS depends on 
the size, shape of the system, and filling factor v, such that the wavelength 
decreases as v approaches an even integer. The coupling of the states of higher 
Landau bands into the lowest band by the Coulomb interaction of the 2DEG 
is essential for the fine structure of the SIS. 

Even though we have been using a restricted HFA here (total angular 
momentum and spin are good quantum numbers) different results can be 
attained by choosing different initial spin configurations. In Fig. [?] we show 
three stable states with higher energy than the ground state seen in Fig. ^c. 
It is interesting to note that the state with no crossing of spin bands is not 
the ground state. 

The exact shape of the SIS does strongly depend on the confining potential 
and, thus, also the size of the system. As was noted earlier the LB's do not 
twist when [i is crossing higher LB's and the spin splitting is opening up, but 
the uneven opening up produces strong modulation of the spin densities. To 
exclude the possibility that numerical deficiencies are causing the twisting of 
the Landau bands we have tested the stability of the spin-density structures 
by increasing the number of basis states included in the numerical calculation 
and tested different schemes in attaining the convergence of the self-consistent 
problem. No visible changes in the ground state properties were observed. On 
the other hand, the exact shape and formation of the SIS does depend on the 
size of the system emphasizing that we are observing a confined spin-density 
wave (SDW) hereSS 

Two possible problems associated with the HFA come to mind. First, 
the HFA may lead to a ground state that is quite different from the physical 
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one due to the strong exchange force that may be reduced in better approx- 
imations where higher order correlation effects or impurity broadening to a 
high order are included.il It is thus, very reassuring that this type of spin 
inversion and formation of a SDW has been observed in models employing 
the local density approximation (LDA) where the SDW has been observed for 
different approximations of the correlation effects Jl The on-set of the SDW 
is also found to depend on the amount of collision broadening of the LB's, 
but neither the broadening nor the correlation effects prevent it J^l The spatial 
correlation of the 2DEG in two approaching finite-size layers for the common 
filling factor of unity is quite similar to the formation of the SIS here. The 
layer index can be treated as isospin for vanishing separation and the nu- 
merical diagonalization of the many-electron Hamiltonian in a large subspace 
of noninteracting many-electron states includes, in principle, all correlation 
effects in the model to a high degree of accuracy.0 

An important difference of the present SIS in the two-dimensional plane 
to the SDW parallel to B investigated by Brey and Halperini is the fact that 
the wavelength of the present modulation varies strongly with v. This is 
caused by the strong dependence of the effective interaction, or the screen- 
ing, in the 2D plane on zaIU'EHH The SDW found by Brey and Halperin 
has strong reassemblance with the more "traditional one" known in ID elec- 
troninc systems.! The notation SIS is, therefore, used here to emphazise this 
difference. 

The region of filling factors when the electrons are not fully spin polarized 
yet (1 < v < 2) but the system has not entered the regime of the integer 
quantum Hall effect with the lowest LB filled {y = 1) has attracted much 
interest lately. It has been shown that in absence of Zeeman energy the 
lowest energy charged excitations at v = 1 are skyrmions, spin textures with 
a unit winding number in two dimensions. 0E1 At large g the quasi-particles, 
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analogous to the single particles, have unit charge ±e and spin half, s = ±1/2, 
but as g is reduced to zero the excitation gap survives and the size of the 
quasi-particles diverges with the spin becoming macroscopic - skyrmions.il 
This effect has also been studied in double-layered electron systems when 
the distance between the layers, each having no spin degree of freedom, is 
reduced since these models can be mapped directly onto the spin system 
identifying the layer index as an isospin.00 It has also been found that these 
spin textures might eventually dominate the ground state properties at filling 
factors 1 < v < 2. The SIS's found in the present work are not related 

to the skyrmions observed in such regime of filling factors, but the skyrmions 
and the SIS's may coexist, which emphasizes the very complex and interesting 
structure of the 2DEG in such a regime. 

The spin-density modulation was found to cause clear signs in the FIR- 
absoption of the confined 2DEG. The signs may be weak since they partly 
reflect single-electron transitions rather than collective oscillations and they 
may be in the low frequency part of the spectrum most difficult to measure, 
but the final word about the appropriatness of the HFA or the LDA for the 
current model will come from experiments. 
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FIGURES 

FIG. 1. The ground state HFA energy spectra and chemical potential \x (horizontal line) for 
(a) N s = 82, (b) N s = 62, (c) N s = 52, and (d) N s = 48 electrons in the system. T = 10.0 K, 
R = 1000 A, Uo = 60 meV, and B = 3.0 T. Crosses represent the s = +| electrons and diamonds 
represent s = — ^. GaAs bulk parameters: m* = 0.067m-o, k = 12 A, g* = —0.44. 

FIG. 2. The ground state energy spectra and chemical potential \i (horizontal line) for (a) 
N s = 42, (b) N s = 38, (c) N s = 34, (d) N s = 30, (e) N s = 26, and (f) N s = 22 electrons in the 
system. T = 4.0 K, R = 1000 A, Uo = 60 meV, and B = 3.0 T. Crosses represent the s = +^ 
electrons and diamonds represent s = — ^. 

FIG. 3. The ground state electron density n s {r) for s = —1/2 (solid) and s = +1/2 (dashed) 
in the case of T = 4.0 K. Other parameters are as in Fig. ||. 

FIG. 4. HFA ground state energy spectra and chemical potential for 30 electrons in a 
parabolic confinement potential at T = K for (a) hu c = 3.0 meV, (b)huj c = 5.0 meV, and 
(c) huj c = 7.0 meV. Only one LB is considered. Confinement frequency huio = 5.0 meV. Crosses 
represent the s = +5 electrons and diamonds represent s = — \. Other parameters are as in Fig. 0. 

FIG. 5. HFA ground state energy spectra and chemical potential /U for 30 electrons in a 
parabolic confinement potential at T = 1.0 K. 1 LB is used for the calculation of the top sub- 
figure, 2 for the center one, and 3 for the bottem one. hu; c = 8.0 meV, and the confinement 
frequency Tiujq = 4.0 meV. Crosses represent the s = +| electrons and diamonds represent s = — ^. 
Other parameters are as in Fig. [p. 

FIG. 6. The FIR-absorption P{E) vs. E/E c {E c = huj c ) and the number of electrons iV s for 
N p = ±1 for HA (left), HFA (center), and for N p = -1 in the HFA (right). T = 4.0 K and other 
parameters are as in Fig. |J 

FIG. 7. Energy spectra and chemical potential fi (horizontal line) of several stable excited 
states for N s = 34. Same parameters as in Fig. |2[ 
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